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Abstract 

In this Letter we research the space charge hmiting current value at which the 
oscihating virtual cathode is formed in the relativistic electron beam as a function 
of the external magnetic field guiding the beam electrons. It is shown that the space 
charge limiting (critical) current decreases with growth of the external magnetic 
field, and that there is an optimal induction value of the magnetic field at which the 
critical current for the onset of virtual cathode oscillations in the electron beam is 
minimum. For the strong external magnetic field the space charge limiting current 
corresponds to the analytical relation derived under the assumption that the motion 
of the electron beam is one-dimensional [High Power Microwave Sources. Artech 
House Microwave Library, 1987. Chapter 13]. Such behavior is explained by the 
characteristic features of the dynamics of electron space charge in the longitudinal 
and radial directions in the drift space at the different external magnetic fields. 
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1 Introduction 

Beams of charged particles have great importance for understanding physi- 
cal properties of plasmas, as well as for their technological applications. High 
power relativistic electron beams are used in plasma heating, inertial fusion, 
high-power microwave generation, and other. One of the most important direc- 
tions of research in the high-power vacuum and plasma electronics is the study 
of the virtual cathode oscillators (vircators). The vircator is a microwave gener- 
ator whose operation is based on oscillations of the virtual cathode ( VC) in the 
electron beam with the overcritical current [1-3]. In recent years vircators have 
bean widely studied experimentally [3-10], by computer simulation [3,11-17] 
and analytical analysis [18-21]. Microwave radiation in vircators is generated 
when the current of a charged particle beam exceeds the critical (so-called 
space charge limiting) current Iscl [1)2]. Therefore, in order to understand 
the physical processes occurring in beams with the VC, it is important to 
determine the critical current Iscl for the formation of the oscillating VC in 
the electron beam. 

The problem of the limiting current for a beam injected into a diode gap was 
first treated in works of Child and Langmuir [22-24]. This problem was later 
considered by Pierce [25], who considered one- dimensional dynamics of the 
electron beam in the planar diode gap with a fixed neutralizing background 
(the so-called Pierce diode) and stated that the increase of the beam current 
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in Pierce diode can lead to instability (which came to be called the Pierce 
instability [26-29]) with the resulting formation of the VC [1]. Further studies 
presented that the development of the Pierce instability in the diode gap is 
accompanied by the various nonlinear oscillatory phenomena, such as chaotic 
spatiotemporal oscillations and the formation of coherent structures (see, for 
example, [26,30-32]). Let us note, that the studies by Child and Langmuir 
were later generalized numerically to beams and drift spaces having a more 
complicated geometry [33]. 

Detailed analysis of the critical current values of relativistic electron beams 
with different geometries was given firstly by Bogdankevich and Rukhadze 
[18]. The work of Genoni and Proctor [19] produces later a more accurate 
relation for space charge limiting current (see also [1]). These results were 
obtained under the assumption that the motion of the electron beam is one- 
dimcnsional (or similarly, that the guiding electrons external magnetic field 
is infinitely strong). The theoretical condition for the space charge limiting 
current assuming one-dimensional electron motion is [1,2,18,19] 

Isc, = ^ ' 70 = - (1) 

e l + 2ln{R/ny ^{l-{v,/cf)' ^ ^ 

where R is the drift tube radius, is the radius of solid beam, 70 is the value 
of relativistic factor at injection in the drift space, vq is the electron velocity 
at injection in drift space, c is the light speed, e is the electron charge, and m 
is the electron rest mass. 

Nevertheless, the problem of how the critical current for the onset of the 
VC in the electron beam depends on the external axial magnetic field has 
been not completely studied. It appears therefore worthwhile to examine this 
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issue in order to gain clearer and deeper insight into the physical processes 
occurring in charged particle beams with the VC. It should be noted that in 
Refs. [14, 17,34] the complex nonlinear dynamics of the electron beams with 
the VC in the finite external magnetic field was studied, but in those papers 
the problem of the dependence of the critical current for the onset of the VC 
oscillations on the parameters of the system, was not analyzed. 

In the present Letter we report the results of numerical studies of how the 
critical current Iscl for the formation of the unsteady oscillating VC in the 
weakly relativistic electron beam depends on the external axial magnetic field. 

The Letter is organized as follows. In Section [2] the mathematical model of 
studied electron-beam system is briefiy discussed. In Section [3] the results of 
investigations of how the external magnetic field affects the critical current at 
which the VC oscillations arises in the electron beam are presented. In Section 
4 we discuss the physical processes that accompany the formation of the VC in 
the electron beam and govern the dependence Iscl{B) described in Section [3] 
for different values of the external magnetic field. In conclusion, we summarize 
the main results discussed in our Letter. 

2 General Formalism 

The simulation model will be briefiy described in this section. The drift space 
for the electron beam is the closed finite-length cylindrical waveguide region 
of length L and radius R with grid electrodes at both ends. An axially- 
symmetrical monoenergetic solid electron beam with the current I and radius 
Rb is injected at the velocity Vq through the left (entrance) grid and than 
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can be is extracted through the right (exit) grid or escapes to the side wall 
of the drift tube. An external uniform guiding electrons magnetic field with 
induction B is applied along the waveguide axis. 

We consider a time-dependent 2.5D model in which the dynamics of the elec- 
tron beam in the drift space is described by solving a self-consistent set of 
Vlasov and Poisson equations [35,36]. The Vlasov kinetic equation for electron 
beam motion analysis is solved numerically by the particle method [35,37,38]. 
The potential distribution in the drift waveguide can thus be easily obtained 
by numerical solving the Poisson's equation in cylindrical geometry [35]. 

The model used in our study does not provide a correct analysis of the problem 
of the critical current of electron beams with the normalized velocity of (3q — 
Vq/c > 0.5, because it does not take into account the self-magnetic field of the 
electron beam and the associated beam pinching effects. Therefore we restrict 
ourselves to studying the critical currents of weakly relativistic electron beams. 

The equations describing the dynamics of the weakly relativistic electron beam 
are formulated in terms of the following dimensionless variables of the potential 
0, the space charge field E, the electron density p, the electron velocity v, the 
spatial coordinates z and r, and the time t: 

^' = (^o/^o) E' = (v^/Lrio) E, B' = (vo/Lrjo) B, p' = poP, 

(2) 

v' — VqV, z' — Lz, r' — Lr, t' — (L/vq) t, 

Here the prime means the dimensionless values, 770 is the specific charge of 
electron, vq and po a-re the velocity and space charge density of the electron 
beam at the entrance to the system, respectively, and L is the length of the 
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interaction space. 



As mentioned above, unsteady processes in the electron beam injected into the 
drift space were simulated numerically by the particle method. In cylindrical 
geometry particles have the form of a charged ring. For each of those charged 
particles the equations of motion were solved with allowance for the weakly 
relativistic velocity of the beam electrons. In terms of dimensionless values 
([2]), the equations of motion in cylindrical coordinates are written as 



dt'^ dt dt 



d'^Zi 



r]{zi,9i,ri)E^ 



dt ' 
1 = 1,... No, 



(3) 

(4) 
(5) 



where 



r]{zi,6i,ri 



PI 



' drj 
dt 



+ Ti 



de^ 

dt 



+ 



' dzj 
dt 



(6) 



Here, Zi, r-i, and di are the longitudinal, radial, and azimuthal coordinates of 
the charged particles, respectively, Ez and are the longitudinal and radial 
electric field components, respectively, B = B^ is the longitudinal component 
of the external magnetic field, whose radial component is assumed to be zero 
Br = 0, /3o = fo/c, where fo is the electron beam velocity at the entrance to the 
system and c is the light speed. The fields do not have azimuthal components 
because the system is axially-symmetric. The subscript i denotes the number 
of particles and Nq is the full number of particles used to model the charged 
particles beam. 



The potential distribution in the interaction space was calculated self-consistently 
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from Poisson's equation 



1 dip ^ (fip ^ dip^ 2 
r dr dr"^ dz^ ' 



where 



a 



/ I I \ 1/2 



is the dimensionless control parameter which depends on the beam current as 
a ~ \/l and is proportional to the length of the drift tube space as a ~ L. 
Poisson's equation ([7]) was solved with the boundary conditions 



ip(z = 0, r) = 0, ip{z = 1, r) = 0, (p{z, r = R) = 0, (9) 
= 0. (10) 



dip 
dr 



r=0 



Condition implies that drift region bounded by the conducting cylindrical 
surface of the radius R held at a zero potential. Condition ffTOj) means that 
the radial electric field at the symmetry axis r = is zero due to the axial 
symmetry of the interaction space. 

For each particle, equations of motion ©-([S]) were integrated numerically by 
the second-order leap-frog method [35]. At each time step Poisson's equation 
([7]) was solved on the two-dimensional mesh in cylindrical coordinates. In order 
to reduce the mesh noise, the space charge density on the mesh was calculated 
using a bi-linear weighing procedure for particles (the particle-in-cell (PIC) 
method) [35]. 
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3 Critical current dependence on the external magnetic field 

Let us consider the results of numerical study of how the external magnetic 
field affects the space charge limiting current Iscl at which the VC oscilla- 
tions arises. The value of the critical current is determined by the calculation 
of the trajectory of the charged particles in the phase space (r, v), where r is 
a radius-vector and v is a velocity of the charged particle. The moment when 
virtual cathode appears corresponds to reflection of a part of the electron 
beam backwards to the injection plane z = 0. In this case, appearance of the 
charged particles with negative longitudinal components of velocity {vz < 0) 
takes place. Onset of the oscillations of the virtual cathode is investigated 
by the analysis of the dynamics of the space charge electric field (both lon- 
gitudinal E^, and radial Er components) in the drift space. The moment of 
appearance of charged particles refiections in the beam (virtual cathode for- 
mation) corresponds to the moment of onset of the electric field oscillations in 
the system. Let us note that if an injected current exceeds the critical current 
value insignificantly then the duration of the time interval required for virtual 
cathode formation and the refiected particle occurrence grows considerably 
(see [5]). Therefore we took into account this fact in the numerical simulation 
by means of choice of the duration of the numerical simulation large enough 
for the beam current being close to the value of the critical current. So, in 
this case the duration of the simulation exceeded 200Tp, where Tp — 27r/a;p 
is the typical time scale, cUp is the plasma frequency of the electron beam. 
This duration corresponds approximately 100 periods of the virtual cathode 
oscillations. 

We begin by analyzing the case in which the ratio of the initial beam radius 
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rfc at the entrance to the interaction space to the radius R of the interaction 
space itself a = rt/R is fixed. We call a the filling parameter, i.e. the param- 
eter that characterizes the extent to which the drift space is filled with the 
electron beam. If it is not stipulated especially, the ratio of the drift tube ra- 
dius to its length is assumed to be constant and is equal to R/L = 0.25. Fig. [1] 
shows how the normalized critical beam current Iscl depends on the external 
magnetic field B for a = 0.5. The critical beam current is normalized to its 
magnitude Iscl at a zero external magnetic field (-8 = 0). In Fig. [1] different 
curves correspond to different initial velocities = vq/c of the electron beam 
injected into the drift tube. We can see that the dependence Iscl{B) behaves 
differently in two characteristic ranges of the magnetic field B. At weak mag- 
netic fields the critical beam current decreases monotonically as the external 
magnetic field B increases. At strong external magnetic fields the critical cur- 
rent of the electron beam increases monotonically with magnetic induction. 
Consequently, there is a certain magnetic field -Bmin at which the critical beam 



Fig. 1. Normalized critical (space charge limiting) current Iscl/ Iscl{B = 0) vs 
external magnetic field B for o" = 0.5 and for the different electron beam velocities: 
(1) /?o = 0.1, (2) 0.3, and (3) 0.5. In the frame the optimum external magnetic field 
-Bmin vs normalized electron velocity /?o is shown. 
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Fig. 2. Critical (space charge limiting) current IscL [kA] vs normalized electron 
beam velocity (3q for a = 0.5 and strong magnetic field B = 2.5 kG (points o). 
The solid line corresponds to the analytical relation ([T]) of the space charge limiting 
current. 

current is minimum. This result is quite important because it provides the way 
to determine the optimum magnetic field for the onset of microwave genera- 
tion in the virtual cathode oscillator at the minimum current of relativistic 
electron beam. Note, that at strong magnetic fields the critical beam current 
increases and asymptotically approaches to the analytical relation ([1]) of space 
charge limiting current which was obtained for an infinitely strong external 
magnetic field (and, as consequence, one-dimensional motion of the electron 
flow) [1] . For verification of the results of our simulation in Fig. [2] the analyt- 
ical (Eq. ([T])) and numerical dependencies of the critical current Iscl on the 
normalized electron beam velocity /?o for strong magnetic field B = 2.5 kG are 
shown. Comparison of the presented dependencies shows good agreement of 
the known analytical results and our numerical calculations. 

The dependencies shown in Fig. [1] were calculated for different values of the 
electron beam velocity (3q at the entrance to the drift tube. We can see from 
Fig. [H that for the non-relativistic and relativistic electron beams under con- 
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Fig. 3. (a) Normalized critical electron beam current Iscl/Iscl{B = 0) vs external 
magnetic field for different values of the filling parameter a (the extent to which the 
drift tube is filled with the beam). Curve 1 corresponds to cr = 0.2, 2 - 0.5, and 3 - 
0.9. (b) Normalized critical current vs. filling parameter a for different inductions of 
the external magnetic field: (1) B = 0, (2) B = B^i^ = 0.6 kG, and (3) B = 2.5 kG. 
The electron beam velocity at the entrance to the drift tube is such that /?o = 0.5 

sideration here the dependence of the critical beam current Iscl on the ex- 
ternal magnetic field B is qualitatively the same for different beam velocities. 
For any velocity of the injected beam the critical current decreases with in- 
creasing B in the range of the weak external magnetic fields, and, otherwise, 
it increases monotonically in the range of the strong magnetic fields. However, 
the optimal magnetic field -Bmin depends on the electron beam velocity Pq. 
The dependence of -Bmin on (3q is shown in the frame in Fig. [T] from which we 
can see that the optimum external magnetic field increases linearly with the 
velocity /5o of the injected electron beam. 

Let us now consider the dependence of the critical current for the formation 
of the VC oscillation on the external magnetic field for different values of the 
filling parameter a and the ratio R/L of the drift tube radius to its length. 
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Fig. 4. (a) Normalized critical electron beam current Iscl/ Iscl{B = 0) vs external 
magnetic field for different values of the ratio of the drift tube radius to its length for 
fixed beam radius r^/L = 0.25. Curve 1 corresponds to R/ L = 1.3, 2 - 1.0, 3 - 0.5, 
4 - 0.4, and 5 - 0.3. (b) Normalized critical current vs. the ratio R/ L for different 
inductions of the external magnetic field: (J) B = 0, {2) B = i?min = 0.6 kG, and 
(3) B = 2.5 kG. The electron beam velocity at the entrance to the drift tube is such 
that Po = 0.5 

Fig. [3](a) presents the corresponding dependencies Iscl{B) for different a 
values and for the fixed electron beam velocity I3q = 0.5 at the entrance to 
the drift tube. We can see, that the current required for the formation of the 
oscillating VC increases with the filling parameter. In Fig. [3t^b) the critical 
beam current dependencies on the filling parameter a for different values of 
the external magnetic field are shown. It can be seen that the critical beam 
current Iscl increases monotonically with a, and that for different external 
magnetic fields the dependence Iscl{o') remains essentially the same. 

Fig. 111(a) shows the dependencies Iscl{B) for different values the ratio R/L 
of the drift tube radius to its length and for the fixed electron beam velocity 
(3q = 0.5. In this case the current required for the formation of the oscillating 
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VC decreases with the growth of the R/L value. In Fig.[3](b) the critical current 
dependencies on the geometry parameter a for different values of the external 
magnetic field are shown. It can be seen that the critical beam current Iscl 
decreases monotonically with R/L, and that for different external magnetic 
fields the dependence Iscl{R) remains essentially the same. 

It should be noted that the optimal magnetic field -Bmin, at which the critical 
beam current is minimal, does not essentially depend on the filling parameter 
a and geometry parameter R/L and is determined only by the velocity (3q. 
For different values of the parameters a and R/L the critical beam current 
accepts minimal value for the optimal external magnetic field B = -Bmin (see 
curves 2 in Figs. [3](b) andlD^b)). 

4 Physical processes in the electron beam with overcritical current 
at the external magnetic field 

Let us consider the physical processes that accompany the oscillating VC 
formation in the electron beam and govern the dependence Iscl{B) described 
in previous Section [3] for different inductions B of the guiding magnetic field. 

We begin by analyzing the dependence of the beam current distribution in 
the system on the external magnetic field. To do this, we consider how the 
numbers of large charged particles that escape from the drift space through its 
entrance and exit planes and through its side surface depending on the external 
magnetic field. Fig. [5] shows how the normalized numbers N of particles that 
are reflected from the VC and escape from the drift space through the entrance 
grid (curve 2 in Fig. [5]) and through the side surface (curve 1) and that pass 
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through the VC and escape from the drift space through the exit grid (curve 
3) depend on the axial magnetic field in the drift space. The calculations 
were carried out for the filling parameter cr = 0.5 and for the beam velocity 
Po = 0.5. 

From Fig. O we can see that the escape dynamics of the charged particles 
from the drift space differs radically in the cases of the weak and the strong 
values of the external magnetic field. In the case of the weak magnetic field 
B < 0.5 kG, a large number of charged particles escape from the drift space 
through its side wall. The number of such particles, that are reflected from 
the VC and escape from the system in the radial (transverse) direction, is 
maximum at zero external magnetic field and decreases monotonically with 
magnetic induction increasing (see curve 1 in Fig. At the same time, the 
number of particles that are reflected from the VC in the longitudinal direction 
and escape through the injection plane z = 0, as well as the number of particles 
that pass through the oscillating VC and escape from the drift space through 
the exit grid, increases with the external magnetic field (see Fig. [5], curves 2 
and 3 respectively). In Fig. [T] such behavior of the number of charged particles 
that are reflected from the VC in the longitudinal and transverse directions 
corresponds to the range B < -Bmin in which the critical beam current decreases 
with increasing external magnetic field (cf. the curves in the range B < 0.5 kG 
in Fig. [T] versus the curves in the same range in Fig. 

For stronger magnetic fields, B > 0.5 kG there are almost no charged particles 
that are reflected from the VC in the transverse direction and escape radially 
from the drift space through its side wall (see Fig. [5]). The number of particles 
moving in the longitudinal direction and escaping from the drift space through 
the entrance and exit grids remains unchanged with the increasing of the 
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Fig. 5. (a) Normalized number N/Nq of charged particles escaping from the drift 
space through the side wall (curve 1) and through the entrance and exit grid elec- 
trodes (curves 2 and 3, respecticely) versus external magnetic field for Pq = 0.5 and 
a = 0.5. The value Nq is the total number of injected charged particles 

external magnetic field (see Fig. [5], curves 2 and 3). A comparison between 
Figs. [Hand [5] shows that the magnetic field at which the transverse dynamics of 
the charged particles reflected from the VC becomes unimportant (i.e. ^ 0) 
and is close to the optimal magnetic field -Bmin at which the critical space 
charge limiting current is minimum. 

Such behavior of the electron beam allows us to suggest that at the weak 
external magnetic field the formation and dynamics of the VC oscillations 
differ from those at the strong magnetic field. In the weak magnetic field 
B < Brain, which poorly confines the electrons, the transverse dynamics of 
the electron beam plays a predominant role in the drift space. For the strong 
external magnetic field that the current through the side surface bounding 
the drift tube in the radial direction is zero the longitudinal dynamics of the 
space charge of the beam is dominated in the analyzed system. Let us note that 
cyclotron rotation of electrons in the magnetic field is observed simultaneously. 
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Fig. 6. Space charge density distributions in the VC region, averaged over the char- 
acteristic period of VC oscillations for the filling parameter a = 0.5 and different 
inductions of the external magnetic field: (a) B = 0, (h) B = i?min = 0.6 kG, and 
(c) B = 2.0 kG. The arrows show the initial radial beam boundary. Longitudinal 
coordinate is shown on the abscissa and radial coordinate is shown on the ordinate. 
The gray color intensity is proportional to the absolute value of the space charge. 
The scale from the right side of the figure shows the values of the space charge 
density. 

At the intermediate magnetic fields B ~ Smin the behavior of the system is 
governed by both the transverse and longitudinal dynamics of the space charge 
in the VC region. 

Let us examine in more detail the behavior of the electron beam in the VC 
area at different inductions of the external magnetic field. The space charge 
density distributions in the area of VC oscillations {z e (0,0.2)) averaged 
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over the characteristic period of VC oscillations for different inductions B is 
shown in Figs. [6l In these figure the absolute values of the dimensionless space 
charge density \p{r, z) \ at different points (r, z) in the drift space are shown by 
shades of gray. The maximum in the space charge density (the dark region in 
the distributions \p{r,z)\) corresponds to the VC region, where the electrons 
are stopped and start moving back toward the injection plane or toward the 
side wall of the drift tube. The electron beam is decelerated by the space 
charge of the VC and the electron velocities within it are close to zero. So, the 
space charge density \p\ achieves the maximum value in the VC region. 

From the distributions shown in Figs. [6]^a), which were calculated for the zero 
external magnetic field, we can see that in the VC region the electron beam 
greatly expands in the transverse (radial) direction, so the beam space charge 
density near the injection plane is nonzero over the entire cross transverse 
section of the drift tube. For the weak external magnetic field the space charge 
density distribution is qualitatively the same. This behavior of the electron 
beam can be explained as follows. A dense electron bunch formed in the VC 
region (and, accordingly, decelerates the electron beam) gives rise to fairly 
strong space charge forces (the Coulomb repulsion between the electrons), 
under the effect of which the electrons acquire the transverse velocity, rapidly 
escape in the radial direction, and are lost at the side cylindrical wall bounding 
the drift tube. 

Such a transverse electron dynamics under the effect of space charge forces in 
the absence of the external magnetic field or in the weak external magnetic 
field, which does not restrain these forces, results in substantial radial expan- 
sion of the beam. As a consequence, in the weak external magnetic field most 
of the beam electrons in the VC region move in the transverse (radial) direc- 
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tion with respect to the drift tube axis. The radial beam expansion reduces the 
space charge density in the VC and, accordingly, increases the space charge 
limiting current Iscl for achieving the space charge density required for the 
formation of the oscillating VC. In this case some of the electrons begin to be 
reflected to the entrance grid and are turned to the side surface of the drift 
tube. 

As can be seen from Fig. [6](a), the electron bunch in the VC region extends 
over the entire transverse cross section of the drift tube. In this case, the 
space charge density is maximum (|p| > 0.25 in dimensionless units) within 
the radial region whose radius exceeds the initial beam radius (shown by 
the arrow in Fig. [n](b)). The radial region of maximum space charge density 
is displaced in the propagation direction of the injected electron beam. This 
indicates that the distribution p(r, z) in the VC region is nonuniform in both 
the radial and longitudinal directions. 

For the small filling parameter a the system behaves in the similar manner. In 
this case the lower critical current of electron beam is determined by the higher 
space charge density of the electron beam because of the smaller beam cross- 
sectional area for the fixed parameter a which is proportional to the injected 
beam current /. Consequently, the space charge density level required for the 
formation of the VC oscillations in the beam of larger radius is achieved when 
the electron beam current is higher than that of the beam with smaller radius. 
This conclusion, which is confirmed by numerical simulations (see Fig. [3]), is 
valid for an arbitrary induction B of the external magnetic field. 

Let us note, that for the weak magnetic field the electron bunch in the VC 
region oscillates mainly in the radial direction. In this case the maximum space 
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charge density and the spatial position of the VC vary only slightly with time, 
in contrast to the case of strong external magnetic fields, in which the electron 
dynamics is essentially one-dimensional (see, e.g., [1,16,39]). Radial oscillations 
can occur with the lower beam space charge densities than those at which the 
VC can execute longitudinal oscillations. 

A stronger magnetic field B prevents the beam electrons moving in the radial 
direction and imparts the azimuthal velocity to them. As a result, the number 
of electrons escaping from the drift space radially through the side wall of the 
system decreases as the magnetic induction B increases (see Fig. [5]). Fig. [6](b) 
shows the averaging space charge density distributions in the VC region for 
the optimal magnetic field B = -Bmin- It is possible to see that the external 
magnetic field restricts radial electron motion, increasing the space charge 
density in the VC region (cf. the gray scale on the right of Figs. |6]^a) and[6]^b)). 
Consequently, as follows from Fig. [1] decreasing the critical current at which 
the electrons begin to be reflected and the oscillating VC arises in the beam. 

Fig. El(c) shows that in the strong magnetic field such that the beam dynamics 
close to one-dimensional {B ^ -Bmin), the VC (the region of maximum space 
charge density) ranges almost entirely within the initial beam radius. The 
initial beam radius is indicated by the arrow in Fig. O In this case radial 
oscillations of the space charge density are not observed. At strong magnetic 
fields the oscillatory behavior of the VC agrees well with the VC oscillations 
in the one- dimensional case described in detail in Refs. [21, 28, 31] for the 
ID model of the electron beam with the VC. As consequence, VC oscillations 
behavior and the space charge limiting current weakly depends on the external 
magnetic field in the range B ^ i?min- 
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Based on the fact that the space charge dynamics differs radically in the cases 
of the weak and the strong external magnetic field, we may speak of two 
different types of oscillating VC. In the weak magnetic field, the formation of 
the VC and its oscillations are governed primarily by the transverse dynamics 
of the beam electrons under the effect of the repulsive space-charge forces, 
which substantially expands the beam in the radial direction. In this case we 
may speak of the transverse VC in which the space charge density oscillates 
predominantly in the radial direction. The situation is radically different for 
the strong external magnetic fields when the longitudinal electron dynamics 
is dominated and the space charge density oscillates predominantly in the 
longitudinal direction. In this case we may speak of the longitudinal VC in 
an electron beam with an overcritical current. The behavior of such beams 
with the longitudinal VC oscillations has been studied in more detail by the 
numerical simulation [11-13,15,16]. Such oscillating regimes are typical of 
vircators operating with a strong external magnetic field [3, 15]. Accordingly, 
the critical beam current is described by the familiar relationship ([1]) [1]. 

At intermediate external magnetic fields B ~ i?min competition between the 
longitudinal and transverse oscillations of the space charge density in the 
electron beam with the VC gives rise to the minimum in the dependence 
Iscl{B) of the space charge limit current on the external magnetic field. 

5 Conclusion 

Thus, the influence of the external magnetic field on the space charge limiting 
current corresponding to the onset of the oscillating VC in the electron beam 
is reported for the first time in this Letter. It is shown that the critical current 
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decreases when external magnetic field increases. There is the optimal mag- 
netic field at which the space charge limiting current is minimum. For stronger 
external magnetic fields the critical current is described by the familiar rela- 
tion ([1]) [1] obtained under the assumption that the motion of the electrons 
is one-dimensional. Study of the physical processes in the electron beam with 
the virtual cathode presents that the dependence of the critical current on 
the external magnetic field is governed by the characteristic features of the 
longitudinal and radial dynamics of the electron space charge in the drift tube 
at different inductions of the axial magnetic field. 

In the weak external magnetic field the major role is played by the transverse 
beam dynamics (i.e. the two-dimensional effects associated with the Coulomb 
repulsion of the electrons in the radial direction). As a consequence, the space 
charge density in the VC area reduces, so the oscillating VC can arise in 
the electron beam with the higher current in comparison with that at strong 
magnetic fields. When the external magnetic field is strong the transverse 
beam dynamics is less important and the behavior of the VC is governed 
primarily by the longitudinal oscillations in the electron beam. 
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